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Abstract 

We investigate the low-temperature behavior of ferromagnets with a spon- 
taneously broken symmetry 0(3) — ► 0(2). The analysis is performed within 
the perspective of nonrelativistic effective Lagrangians, where the dynam- 
ics of the system is formulated in terms of Goldstone bosons. Unlike in a 
Lorentz-invariant framework (chiral perturbation theory), where loop graphs 
are suppressed by two powers of momentum, loops involving ferromagnetic 
spin waves are suppressed by three momentum powers. The leading coeffi- 
cients of the low-temperature expansion for the partition function are cal- 
culated up to order p 10 . In agreement with Dyson's pioneering microscopic 
analysis of the cubic ferromagnet, we find that, in the spontaneous magneti- 
zation, the magnon-magnon interaction starts manifesting itself only at order 
T 4 . The striking difference with respect to the low-temperature properties of 
the 0(3) antiferromagnet is discussed from a unified point of view, relying on 
the effective Lagrangian technique. 
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I. INTRODUCTION 



In the following presentation, our interest is devoted to the low-temperature behavior of 
ferromagnets, which exhibit a spontaneously broken internal symmetry 0(3) — > 0(2). This 
system has been widely studied in condensed matter physics and it is not our intention to 
contribute to its detailed microscopic understanding. Rather, we want to analyze several 
low-energy phenomena from a unified point of view, relying on the method of nonrelativistic 
effective Lagrangians. We try to understand how the symmetries, inherent in the underlying 
theory, manifest themselves in the partition function at low temperatures. The complex 
microscopic description of the system is taken into account only through a phenomenological 
parametrization, which, in the effective Lagrangian, emerges in the form of a few coupling 
constants. 

Nevertheless, let us first consider the Heisenberg model, which describes the ferromagnet 
on a microscopic level. There, the exchange Hamiltonian T~C , 

H = -J^2S m -S n , J = const, (1.1) 

n.n. 

— * 

formulates the dynamics in terms of spin operators S m , attached to lattice sites m. Note 
that the summation only extends over nearest neighbors and that the isotropic interaction 
is assumed to be identical for any two adjacent lattice sites. For positive values of the 
exchange integral J, the above expression leads to an adequate low-energy description of 
the ferromagnet. 

The interaction between a constant magnetic field H = (0,0, H), H > 0, and the spin 
degrees of freedom is taken into account through the Zeeman term. In the corresponding 
extension of the Heisenberg model,0 

H = Ho - nY,Sn-H, (1.2) 

n 

the magnetic field is coupled to the vector of the total spin. 

Within this model, Dyson evaluated the low-temperature expansion for the partition 
function of a ferromagnet for all three types of cubic lattices.! In particular, for the most 
prominent order parameter of the ferromagnet, the spontaneous magnetization, he obtained 
the following series: 

E(T)/E(0) = 1 - a T§ - a x T^ - a 2 T^ - a 3 T 4 + 0(T*) . (1.3) 

The T 3//2 -term, referred to as Bloch's law, corresponds to free magnons. The next two terms 
involving the coefficients ot\ and «2 are related to the shape of the dispersion curve - or the 
discreteness of the lattice - and describe free magnons as well. Remarkably, the spin-wave 
interaction starts manifesting itself only at order T 4 . 

In order to obtain this result, Dyson had to set up a highly involved mathematical 
machinery. It is our goal to rederive this series within the effective Lagrangian framework 
and to understand Dyson's result by exclusively considering the symmetries of the theory. 
In order to do so, we will first have to establish the momentum power counting scheme for 
this nonrelativistic system, which represents the very basis for a systematic expansion of 



2 



quantities of physical interest. As we will see, it is quite different from the power counting 
scheme in Lorentz-invariant effective theories (chiral perturbation theory). 

Another aspect of the present work concerns the comparison of the low-energy properties 
of 0(3) ferromagnets and 0(3) antiferromagnets. Within the framework of nonrelativistic 
effective Lagrangians, we will, e.g., be able to understand the striking differences in the 
low-temperature expansions for the order parameters: the spontaneous magnetization for a 
ferro magnet and the staggered magnetization for an antiferro magnet, respectively. 

In the effective Lagrangian perspective, the excitations near the ground state, the spin 
waves or magnons in the present case, are interpreted as Goldstone bosons resulting from 
a spontaneously broken internal symmetry. Indeed, the Heisenberg Hamiltonian (EI) is in- 
variant under a simultaneous rotation of the spin variables described by the symmetry group 
G = 0(3), whereas the ground state of a ferromagnet breaks this symmetry spontaneously 
down to H = 0(2): all the spins are aligned in one specific direction, giving rise to a nonzero 
spontaneous magnetization. 

Whenever a physical system exhibits spontaneous symmetry breaking and the corre- 
sponding Goldstone bosons represent the only low-energy excitations without energy gap, 
we do have a very powerful means at our disposal to analyze its low-energy structure: effec- 
tive Lagrangians. The method was originally developed in connection with Lorentz-invariant 
field theories (chiral perturbation theory), admitting a low-energy analysis of the strong 
interaction described by quantum chromodynamics (QCD). The effective Lagrangian tech- 
nique has also proven to be very useful in the investigation of other relativistic systems 
where Goldstone bosons occur, see, e.g., Refs.i~0. The method has been extended to finite 
temperature,^ allowing for a systematic low-temperature analysis of the partition function. 

In condensed matter physics, spontaneous symmetry breaking is a common phenomenon 
and effective field theory methods are widely used in this domain. Only recenth^however, 
has chiral perturbation theory been extended to such nonrelativistic systemsJlHTLa demon- 
strating its applicability to solid state physics as well. The method is based on effective 
Lagrangians which exploit the symmetry properties of the underlying theory, i.e., the Heisen- 
berg model in the present case, and permits a systematic low-energy expansion of quantities 
of physical interest in powers of inverse wavelength. The few applications related to systems 
exhibiting collective magnetic behavior that appeared in the literature so far concern spin- 
wave scattering processes,^ spin-wave mediated non-reciprocal effects in antiferromagnets,^ 
the low-temperature expansion of the staggered magnetization of 0(N) antiferro magnets^ 
and spin waves in canted phases.0 An interesting application concerning SO (5) invariance 
and high-T c superconductivity can be found in RefM. Pedagogical introductions to effective 
Lagrangians are Refs.&M brief outlines of the method can be found in Refs.i. 

As the effective analysis refers to large wavelengths, it does not resolve the microscopic 
structure of a solid and the system hence appears homogeneous. Accordingly, the effective 
Lagrangian is invariant with respect to translations. On the other hand, the effective La- 
grangian is not invariant under rotations, since the lattice structure of a solid singles out 
preferred directions. In the case of a cubic lattice, the anisotropy, however, only shows up at 
higher orders of the derivative expansion:0'0 the discrete symmetries of a cubic lattice thus 
imply space rotation symmetry. In the following, we assume that the ferromagnet exhibits 
this type of lattice structure: the leading order effective Lagrangian is then invariant both 
under translations and under rotations and the corresponding expression for a ferromagnet 



3 



takes the fornMl 



C 2 eff = E a °" + HfiHU 3 - \F 2 d r U l d r W. (1.4) 

In the above notation, the two real components of the magnon field, U a (a = 1,2), have 
been collected in a three-dimensional unit vector U l = (U a ,U 3 ), which transforms with the 
vector representation of the rotation group 0(3). At leading order, the ferromagnet is thus 
characterized by two different low energy constants, £ and F. The first term, which involves 
a time derivative, is related to a topological invariant. Remarkably, due to this contribution 
proportional to the spontaneous magnetization E, the effective Lagrangian of a ferromagnet 
fails to be invariant under the group G = 0(3). Note that such expressions, involving order 
parameters associated with the generators of the spontaneously broken group, would not be 
permitted in Lorentz-invariant effective theories. It represents the main novelty occurring in 
condensed matter physics, where nonrelativistic kinematics is less restrictive than Lorentz 
invariance. 

The associated equation of motion is the Landau-Lifshitz equation, 

d U a + e ajk UU k + je ajk AWU k = , f& = ^H5{ , j = ^ , (1.5) 

which describes the dynamics of ferromagnetic spin waves. Its structure is of the Schrodinger 
type: first order in time, but second order in space. Note that, according to Goldstone's 
theorem, there have to be two real Goldstone fields in the case of a spontaneously broken 
symmetry 0(3) — > 0(2). However, in the nonrelativistic regime, this does not necessarily 
imply that there also have to exist two independent magnon states. Indeed, in the present 
case of a ferromagnet, a complex field is required to describe one particle: there exists only 
one type of spin-wave excitation in a 0(3) ferromagnet, following a quadratic dispersion 
relation, 

u(k) = 7 fc 2 + C(|fc| 4 ). (1.6) 

Accordingly, in the effective description of this nonrelativistic system, two powers of mo- 
mentum correspond to only one power of energy or temperature: k 2 oc u, T. 

The effective Lagrangian method provides us with a simultaneous expansion in powers 
of momentum and of the external field H. The important point is that, to a given order in 
the low-energy expansion, only a finite number of coupling constants and a finite number of 
graphs contribute. Let us now set up the effective Lagrangian for a ferromagnet at higher 
orders in the momentum expansionl! and consider the corresponding power counting scheme. 



II. POWER COUNTING AND EFFECTIVE LAGRANGIAN 

In his pioneering microscopic analysis, Dyson evaluated the temperature expansion for 
the partition function of a cubic ferromagnet up to terms of order T 5 . In our effective 
language, where one power of temperature or energy counts like two powers of momentum, 
this corresponds to an expansion up to order p 10 . We will now show that, to this order in the 
momenta, contributions to the effective Lagrangian involving at most six space derivatives 
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enter and the perturbative expansion requires the evaluation of graphs containing at most 
two loops. 

First of all, note that there are no contributions to the effective Lagrangian with an odd 
number of space derivatives: parity excludes terms like 

c abc e rst d r U a d s U b d t U c , (2.1) 

which involve the antisymmetric tensor e rst . 

Before constructing the relevant pieces C\fp ■ ■ ■ step by step, we have to point 
out an important assumption underlying the present analysis: we assume that the effective 
Lagrangian of a ferromagnet is gauge invariant at sub leading orders. In fact, the discrete 
symmetries of a cubic lattice and time-reversal invariance ensure that there is no topological 
contribution in (see Ref.El) . Topological contributions which may show up in or 
beyond do not affect the conclusions of the present work. 

Gauge invariance then implies that the magnetic field H enters the derivative expansion 
through a timelike covariant derivative, 

D IP = doll* + e ijk f U k , fl = fiH6i . (2.2) 

Note that the magnetic field occurs on the same level as the time derivative and is thus 
counted as a quantity of order p 2 . 

The next-to- leading order Lagrangian is of order p 4 . It involves terms with two time 
derivatives, terms with one time and two space derivatives and terms with four space deriva- 
tives. Time derivatives, however, can be eliminated with the equation of motion,@ 

d U a + e ajk f Q U k + 7e ajk AU j U k = 0, 

such that we end up with the following independent terms:!! 

Cl ff = l 1 (d r U i d r U i ) 2 + l 2 (d r U i d s U i ) 2 + l 3 AU i AU i . (2.3) 

We thus have three effective coupling constants, l±, Z 2 and Z 3 , at next-to- leading order. Note 
that all terms involving the magnetic field have been eliminated with the equation of motion. 
In particular, gauge invariance implies that there is no tree graph from contributing to 
the vacuum energy. 

Clearly, the question now arises, to what order n in the effective Lagrangian C 2 Jjj we 
have to go to evaluate the partition function up to accuracy p 10 . In the following, we will 
show that loop corrections involving ferromagnetic magnons are suppressed by three powers 
of momentum and the effective Lagrangian is needed up to £e//- 

To verify this statement, let us consider a scattering process. The effective Lagrangian 
provides us with an expansion for a multi-magnon scattering amplitude in powers of mo- 
mentum. The elastic magnon-magnon scattering amplitude is obtained by expanding £ e // 
up to order (U a U a ) 2 . As shown in Figure la, the leading term corresponds to the tree graph 
involving the four-magnon vertex from L 2 ^ and is of order F 2 k 2 /E 2 (see Ref.0). The first 
correction comes from the tree graph lb with a vertex involving the next-to-leading order 
Lagrangian £e//- However, we also have to consider loop graphs in the effective framework. 
These corrections are in general infinite and need to be renormalized. 
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la lb lc 

FIG. 1. Feynman graphs describing the leading contributions to the magnon-magnon scattering 
amplitude. The numbers attached to the vertices refer to the piece of the effective Lagrangian they 
come from. Vertices associated with the leading term C 2 ^ are denoted by a dot. 

Consider the magnon-magnon scattering graph lc, where the two vertices come from 
C 2 e ff- The form of the loop graph is 




A loop <x -— T — =- , (2.4) 



where the quantities u and k denote generic energy and three momentum, respectively. Each 
four-magnon vertex from C\jf is of order F 2 k 2 /T, 2 , while the two magnon propagators are 

of order — ^k 2 ). By dimensional analysis, the loop integral must then have the form 

Aioop oc ^ \p 1 5 . (2.5) 
Comparing this expression with the leading tree graph amplitude la, 

Atree OC — p , (2.6) 

we find the remarkable result that, in the effective description of the ferromagnet, loops 
are suppressed by three powers of momentum. This is to be contrasted with the Lorentz- 
invariant framework (chiral perturbation theory), where loops are suppressed by two powers 
of momentum in four dimensions. 

Hence, in order to carry out the expansion for the partition function up to accuracy p 10 , 
we need the Lagrangian up to terms of order p 6 : shows up in a one-loop graph at order 
p 9 . Notice that the only relevant term coming from is quadratic in the magnon field. 
Eliminating time derivatives, we end up with 

C 6 eff = dWA 3 ^. (2.7) 

As in £e//5 terms involving the magnetic field do not show up: they have all been eliminated 
with the equation of motion. More generally, assuming £e//> ^eff an< ^ ^eff t° De g au g e 
invariant, there are no tree-graph contributions from these higher order pieces of the effective 
Lagrangian to the vacuum energy - provided that the magnetic field is the only external 
field considered. 
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III. FINITE TEMPERATURE 



In finite temperature field theory, the partition function is represented as a Euclidean 
functional integral ,EHH 

Tr [exp(-H/T)\ = J [dU] exp ( - Jd 4 xC eff ) . (3.1) 

The integration is performed over all field configurations which are periodic in the Euclidean 
time direction U(x, x 4 + 0) — U(x, x 4 ), with (3 = 1/T. The low-temperature expansion of 
the partition function is obtained by considering the fluctuations of the field U around the 
ground state V = (0,0, 1), i.e., by expanding U 3 in powers of U a , U 3 = \/l — U a U a . The 
leading contribution is of order p 2 and contains a term quadratic in U a which describes free 
magnons. In the presence of a magnetic field, they obey the dispersion relationil 

ifa(k) = jk 2 + fiH + 0{\k\ 4 ) . (3.2) 

The remainder of the effective Lagrangian is treated as a perturbation. Evaluating the 
Gaussian integrals in the standard manner, one arrives at a set of Feynman rules which 
differ from the conventional rules of the effective Lagrangian method only in one respect: 
the periodicity condition imposed on the Goldstone field modifies the propagator. At finite 
temperature, the propagator is given by 

oo 

G(x) = A(x, 2:4 + n/3) , (3.3) 

n = — oo 

where A(x) is the Euclidean propagator at zero temperature, 

„ . . r dk 4 d 3 k e ikx-ik4X4 

J (.27TJ 4 jk 2 - %k 4 + fxH 

Note that the above Green function corresponds to the propagation of a single magnon, 
described by the complex field u = U 1 + ill 2 . 

We restrict ourselves to the infinite volume limit and evaluate the free energy density z, 
defined by 

z = —T lim L~ 3 In [Tr exp(— 7i/T)] . (3.5) 

L— »oo 

Temperature thus produces remarkably little change: to obtain the partition function, one 
simply restricts the manifold on which the fields are living to a torus T in Euclidean space. 
The effective Lagrangian remains unaffected - the coupling constants F, E, h, . . . are tem- 
perature independent. 

It is convenient to use dimensional regularization to evaluate the graphs of the effective 
theory. Unlike in a Lorentz-invariant framework, where the physical dimension d is equal to 
four, we regularize the nonrelativistic propagator only in the three spatial components, 

A ( x ) = 77r~v2 (~) 2 ~~~d ex P ( ~ J~ fiHxt) e(x 4 ) , (3.6) 

(2TT) d ^ X4 2 v 47x4 1 
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and set d — 3 at the end of the calculation. 
For the thermal propagator, we then have 

G ( x ) = jJ-^ (t) 2 £ 4r exp ( - -J— - /iHx n ) G(x n ) , (3.7) 

( 27r ) V 7 7 n^ooxi A ^ Xn 

with 

x n = X4 + n/3 . (3.8) 

Using complex notation for the magnon field, the evaluation of Gaussian integrals simplifies, 
since the following two-point functions vanish: 

(0\T{u(x,x 4 )u(y,y 4 )}\0) = (0| T{u\x, x 4 ) u\y, y 4 )}|0) = 0. (3.9) 

Nonzero contributions to the partition function only result from 

(0\T{u(x,x,)u\y,y,)}\0) = ^A(x-y). (3.10) 

The calculation simplifies further due to space rotation symmetry of the propagator, imply- 
ing that, at the origin, single space derivatives vanish, 

d r G(\x I, x 4 )\ x=0 = 0. (3.11) 

Let us now evaluate the free energy density order by order in the momentum expansion. 

IV. EVALUATION OF THE FREE ENERGY DENSITY 

The relevant Feynman graphs are shown in Figure 2. Depicted are all contributions to 
the free energy density up to order p 10 or, equivalently, order T 5 . 

As we will see, all these contributions to the free energy density exclusively involve prop- 
agators, and space derivatives thereof, to be evaluated at the origin x = 0. It is convenient 
to introduce the following notation: 



x=0 



Gi = [G(x 

[AG(x)] x=o , (4.1) 

where A represents the Laplace operator in three dimensions - no confusion should occur 
with A(x), denoting the zero-temperature propagator. 

The quantity Gi is split into a finite piece, which is temperature dependent, and a 
divergent piece, which is temperature independent, 

G 1 = G\ + G°. (4.2) 

Using dimensional regularization, the explicit expressions are 
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9a 9b 10a 10b 

FIG. 2. Feynman graphs occurring in the low-temperature expansion of the partition function 
for an 0(3) ferromagnet up to order p 10 . The numbers attached to the vertices refer to the piece 
of the effective Lagrangian they come from. Vertices associated with the leading term C? e ^ are 
denoted by a dot. The numbers specifying individual graphs correspond to the power of momentum 
in the derivative expansion. Note that, in the counting scheme for ferromagnetic magnons, two 
powers of momentum correspond to one power of energy or temperature. 



Gr "(27T)-( 7 ) a S 

1 (2ir) d V 



7T\ - 00 e^^HnP 



(n/3) 



d 1 
2 



— exp ( 



X 



4^X4 



e(x 4 ) 



(4.3) 



J x=0 



Likewise, for Ga we have 



G a {2n)d ( 7 ) 2 ( 27 )E 



(n(3) 



(27T) 



d ( 7 ) 



1 



-d 



+ 



X 



X4 



F+ 1 I 27 47 2 x 4 



exp ( — 



.7 



e(x 4 ) 



(4.4) 



x=0 



The temperature- independent pieces G° and G^, as well as propagators involving higher 
orders of space derivatives, are all related to momentum integrals of the form 



d d k (k 2 ^j exp — 7x4 k 2 — x^fiH 
which are proportional to 



m = 0,1,2, 



exp[-x 4 fiH] d 
2— Vim + -) . 



(4.5) 



(4.6) 



In dimensional regularization, however, these expressions vanish altogether: G\, G° A and 
propagators involving more Laplacians do not contribute in the limit d — > 3. Therefore, in 
the evaluation of the above graphs, only the temperature-dependent pieces Gj, G^, . . . are 
relevant. 
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The first contribution to the free energy density, which is of order p 2 , comes from tree 
graph 2. It does not depend on temperature and is given by 



z 2 = —Hfj,H . 



(4.7) 



The order p 5 contribution from one-loop graph 5 is associated with a <i-dimensional 
nonrelativistic free Bose gas. For the temperature-dependent part, we obtain 



1 



8^27 



(4.8) 



n=l 



At order p 7 , the next-to-leading order Lagrangian comes into play. For one-loop 
graph 7, which involves a two-magnon vertex, we get 



2k 



Z7 = 



A 2 G{x) 



x=0 ' 



yielding 



15 h 



167T2 E72 



00 -fj,Hnj3 



(4.9) 



(4.10) 



=1 n2 



The first two- loop graph appears at order p 8 . Remarkably, this contribution is propor- 
tional to single space derivatives of the propagator at the origin and thus vanishes, 



Zg oc 



d r G(x) _ d r G{x) 



x=0 



x=0 



0. 



(4.11) 



At order p 9 , two one-loop graphs show up which involve and £e//> respectively. For 
graph 9a, we get 



Z9a — 



2ci 



A 3 G(x) 



x=0 ' 



leading to the temperature-dependent contribution 

T 105 ci 



'9a — 3 9 

327r2L72 



= 1 «2 



(4.12) 



(4.13) 



Graph 9b is proportional to an integral over the torus T which involves a product of two 
propagators, 



Z9b = 



s 2 



Jd d+1 xA 2 G(x) A 2 G(-x). 



(4.14) 



This expression, however, can be reduced to an x-independent term involving one propagator 
only. To verify this statement, consider the equation for the thermal propagator, 



A_ 7 a + ^G(x) = 5(x), 



(4.15) 
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and take the derivative with respect to the magnetic field. The quantity dG(x) / d(fiH) may 
then be written as a convolution integral over the torus. At the origin, the expression reads 



dG{x) 
d(fiH) 



J x=0 



- Jf +1 y G{-y) G(y) . 



(4.16) 



Inserting Laplace operators at intermediate steps, we obtain the more general result 



^(m+n) 



dG(x) 
d(jjH) 



x=0 



-Jd d+1 y A rn G{-y) A n G(y) 



such that ( |4.14j ) can be written as 

z 9b = 



2Zf 
£ 2 



A 



4 dG(x) 



J x=0 



For the temperature-dependent part of graph 9b, we then get 



945 ll » ^ P -»Hnp 



-9b 



647T2 S 2 7 2 



=1 n2 



(4.17) 



(4.18) 



(4.19) 



Finally, at order p , two-loop graphs with insertions from show up. Graph 10a 



contributes with 



3£ 2 



(8Zi + 6/ 2 + 5/ 3 ) G A G A 



2U 



G 1 A 2 G{-x) 



x=0 



The evaluation of graph 10b amounts to 



ZlOb 



2h 



G x A 2 G{-x) 



x=0 



(4.20) 



(4.21) 



and cancels the second term in zio a . For the temperature-dependent part of order p 10 we 
thus end up with 



3(8/x + 6/ 2 + 5/3) 



T 5 E 



00 ^—fiHnf3 °\ 2 



1287T 3 S 2 7 5 u=1 n2 

Collecting terms, the result for the free energy density up to order p 10 becomes 
1 „s ^ e -» Hnf3 15 h 



(4.22) 



2 2 2 
5^272 



105 



= 1 112 

9Z| 



167T2 £7 



x) —fj,Hnj3 



=1 ^2 



327riS 7 i V 2 7 s 

3(8/x + 6l 2 + 5/3 
128tt 3 S 2 7 5 



T 5 E 



n=l ™ 2 
00 g—fiHnf3 ^ 2 



n=l 



5 

77.2 



(4.23) 



The first term is temperature-independent and originates from tree graph 2. Contributions 
which involve half integer powers of the temperature - T 5 ^ 2 ,T 7 ^ 2 and T 9//2 , respectively - 
arise from one-loop graphs and are all related to the free energy density of noninteracting 
magnons. Remarkably, there is only one term in the above series, the order T 5 contribution 
coming from two-loop graph 10a, which is due to the magnon-magnon interaction. 
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9a 



9b 



0--O- 

10a 10b 10c 

FIG. 3. Feynman graphs occurring in the low-energy expansion of the two-point function for 
an 0(3) ferromagnet up to order p 10 . The numbers specifying individual graphs correspond to the 
power of momentum in the derivative expansion. 

V. INTERLUDE: FREE MAGNONS VERSUS INTERACTING MAGNONS 

Consider the formula for the free energy density of a gas of noninteracting bosons, 

* = z * + j£ys j d " k ln i 1 ~ e ^ Ck)/T ) > ( 5J ) 

where zq is the energy density of the vacuum. With the leading term of the dispersion 
relation, 

u(k) = 7 P+/i# + 0(|£| 4 ), 7 = ^, 

we readily reproduce the dominant one-loop contribution zj oc T 5 / 2 in the free energy 
density of the 0(3) ferromagnet. 

In order to account for subleading terms in the dispersion relation, we have to evalu- 
ate higher order contributions to the two-point function (0| T{u(x, x±) u'(y, 2/4) }|0) of the 
magnon field. The relevant graphs are shown in Figure 3. Depicted are all contributions up 
to order p 10 . Instead of listing individual results for the two-point function, we give the final 
expression for the dispersion relation originating from these graphs: 

u}(k) = ^(F 2 k 2 + E/jlH - 2l 3 k 4 + 2 Cl k 6 ) . (5.2) 

The two new terms correspond to graph 7, graph 9a and graph 9b of Figure 3. Note that 
one-loop graph 8 does not contribute to the dispersion relation. The effect of the three 
remaining one- loop graphs of order p 10 will be discussed in a moment. 
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Inserting the above dispersion relation into formula ( |5.1| ), we reproduce all contribu- 
tions to the free energy density which are associated with one-loop graphs: graphs 5, 7, 9a 
and 9b of Figure 2. These contributions are all related to noninteracting magnons - the 
corresponding coefficients in the dispersion relation are independent of the temperature. 

As far as the evaluation of the three one-loop graphs of order p 10 (Figure 3) is concerned, 
the situation is the following: in the sum of these contributions, only one term remains. It 
originates from graph 10a and leads to an additional term in the dispersion relation 



oc 



8h + 6l 2 + 51. 



G° A k 2 . (5.3) 



However, in dimensional regularization, the quantity G° A , as we have seen, vanishes, such 
that the dispersion relation is not affected by the graphs of order p 10 . 

Nevertheless, we can mimic the magnon-magnon interaction by allowing some of the 
coefficients of the dispersion relation to become temperature dependent. If we replace G° A 
by in (|5.3j ), we obtain an interaction term proportional to T 5 in the free energy density. 
Moreover, we see that the leading order coupling constant 7 is not renormalized at leading 
order (oc Gjk 2 , i.e., order p 8 in the free energy density), but only at next-to-leading order 
(oc G^k 2 , i.e., order p 10 in the free energy density). This observation, implying that the 
two-loop graph 8 does not contribute to the dispersion relation, was already made in Refs.0. 

In this section, we have split the free energy density into a piece corresponding to free 
magnons and a piece representing the magnon-magnon interaction, 

Z £/ree Winter ■ (5-4) 

The free part can be obtained from the free Bose gas formula 

Zfree = Z + J d\ In (l - e^*^) , 

with the modified dispersion relation 

uj(k) = ^F 2 k 2 + - 2l z k 4 + 2 Cl £ 6 j. 

The coupling constants F, S, I3 and c\ are all independent of the temperature. The inter- 
action part originates from two-loop graph 10a which involves a four-magnon vertex from 
C A 
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In order to make the structure of the low-temperature expansion for the free energy 
density of an 0(3) ferromagnet more transparent, we rewrite the series in the form 
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z = z - h T* - hiT* - h 2 Ti - h 3 T 5 + 0(T^), (5.5) 
where the temperature-dependent coefficients hi are given by 
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z = —EfiH , 

I 00 e -^Hnf3 
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1 3 7" / 7 5 
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105 / 9/f \ ~ e-» Hn P 
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In the limit (jiH/T)—*0, these coefficients become temperature independent and the sums 
reduce to Riemann zeta functions, 

^0 — - 3 1 C(|) j 
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VI. SPONTANEOUS MAGNETIZATION 



With the above representation for the free energy density, the low-temperature series 
for some relevant thermodynamic quantities can readily be derived. Since the system is 
homogeneous, the pressure is given by the temperature dependent part of the free energy 
density, 

P = z -z, (6.1) 
and the low-temperature expansion amounts to 

P = h T 5 2 + h x Tl + fr 2 T§ + h 3 T 5 + 0(T^), (6.2) 

with coefficients hi given in ( |5.6|) . 

The first contribution represents the free Bose gas term which originates from one-loop 
graph 5. Subsequent terms exhibiting half integer powers of the temperature are related to 
the shape of the dispersion curve and represent effects due to the discreteness of the cubic 
lattice. Remarkably, there is no T 4 term in this series: two-loop graph 8 does not contribute. 
In the low-temperature expansion for the pressure, the interaction among ferromagnetic spin 
waves only starts manifesting itself at order T 5 . 

The corresponding expressions for the energy density u, for the entropy density s and 
for the heat capacity cy are readily worked out from the thermodynamic relations 
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(6.3) 



Reordering powers of the temperature and taking the limit (fiH/T)— >0, we get 

u = fh T^ + \h x Ti + l~h 2 T°* + 4h 3 T 5 + 0(T%), 
s = fh Tl + \h x T* + l~h 2 Tl + 5/i 3 T 4 + 0(Tf), 
c y = f fco Ti + f /*! rS + f h 2 + 20fc 3 T 4 + 0(T§) , (6.4) 

with coefficients ^ given in ( p.7| ). 

Let us now consider the low-temperature expansion for the spontaneous magnetization, 

With the expression ( |5.5| ) for the free energy density, we obtain the following series: 

E(T)/E = 1 - a Tl - axTl - a 2 T* - a 3 T 4 + 0(T§) . (6.6) 
The coefficients a ra are independent of the temperature and given by 



"0 = ~ 3^ a C( 
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105 / 9Zf \ 
327T2£ 2 72 \ 2 7^ / 
3(8Zi + 6Z2 + 5Z3) f5 . , 3 , 

" 3 = 647r 3E3 7 5 C(a) • ( 6 - 7 ) 

Here comes the appropriate place where we would like to compare our result with Dyson's 
microscopic calculation. In fact, one of the main motivations for his analysis was the follow- 
ing question: at what order in the low-temperature expansion does the spin-wave interaction 
manifest itself in the spontaneous magnetization? Before his rigorous analysis, in which he 
showed that the spin-wave interaction only starts contributing at order T 4 , there appeared 
to be an amazing mess in the literature. At least three different answers were available, all 
of them in contradiction with each other and, as it turned out, also in contradiction with 
Dyson's result: the authors! 4 ! obtained a term proportional to T 2 and the two referencesHH 
ended up with two different contributions of order T 7 / 4 . 

Within the effective Lagrangian framework we clearly confirm Dyson's finding: in the 
spontaneous magnetization, the leading contribution coming from the magnon-magnon in- 
teraction is of order T 4 . Any microscopic information, such as the specification of the type 
of the cubic lattice, is comprised in the numerical values of the effective coefficients a n - they 
involve the coupling constants F, E, Zi, . . . occurring in the derivative expansion of the effec- 
tive Lagrangian C 2 e ^, £ 4 jy and which phenomenologically parametrize the microscopic 
detail. 
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After Dyson's work, various authors tried to derive his results with alternative methods 
within the microscopic framework of the Heisenberg model. In particular, there was emerg- 
ing interest in simplifying his complicated calculation.0 In the course of these investigations, 
however, not all authors could confirm Dyson's findings and, surprisingly, a new term in the 
temperature expansion for the spontaneous magnetization started to haunt the literature: 
various authors^ ended up with an interaction term of order T 3 . However, the T 3 contribu- 
tion turned out to be a spurious effect .0 It emerged due to the approximation methods used 
in these involved microscopic calculations: random phase approximation and decoupling 
approximation in the method of double-time temperature-dependent Green functions. 

We would like to emphasize that there are no such approximations in the effective La- 
grangian framework. The method applies to any system where the Goldstone bosons are the 
only excitations without energy gap. In the effective Lagrangian perspective, the problem is 
approached from a unified and model-independent point of view, based on the symmetries 
inherent in the underlying theory. In the present case of a ferromagnet, it is the sponta- 
neously broken internal symmetry 0(3) — > 0(2) of the Heisenberg model and the nonzero 
value for the spontaneous magnetization, which dictate the structure of the low-temperature 
series considered in this paper. 

Our calculation clearly shows that there is no T 3 term in the low-temperature expansion 
for the spontaneous magnetization of a cubic 0(3) ferromagnet. Using effective language 
this translates into the statement that the contribution from two-loop graph 8 of Figure 2 to 
the free energy density vanishes. As we have pointed out, the reason is due to space rotation 
symmetry of the leading order effective Lagrangian, implying that single space derivatives 
of the propagator, evaluated at the origin, are zero. 

The argument is not restricted to cubic lattices. Consider the leading order Lagrangian 
of the 0(3) ferromagnet: 

e ab d U a U b 
1 + U 3 

Let us replace the last term, which is invariant under space rotations, by the more general 
expression 

- \Fld 1 U i d 1 U i - ^d 2 U%U i - iFldzU'dzU 1 , (6.8) 

in order to account for anisotropics of the crystal lattice. The modified equation for the 
thermal propagator, 

( d 



y () ; 4 - lid( - J2% - 7sdt + nHj G(x) = 5(x) , (6.9) 

again, allows us to eliminate all terms in the evaluation of graph 8 up to contributions which 
involve single space derivatives of the propagator at the origin, 

z 8 = Y, c p[ 9 p G ( x )] x=0 [ d p G ( x )] x=0 > P= 1,2,3. (6.10) 

p XX 

Even though the thermal propagator is no longer space rotation symmetric, it is still invariant 
under parity, such that single space derivatives, evaluated at the origin, vanish. We conclude 
that, also for anisotropic lattices, described by the new contributions ( |6.8| ) in the effective 
Lagrangian, there is no interaction term of order T 3 in the low-temperature series for the 
spontaneous magnetization. 
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VII. ORDER PARAMETERS OF FERRO- AND ANTIFERROMAGNETS: 
SPONTANEOUS VERSUS STAGGERED MAGNETIZATION 



The low-energy behavior of an 0(3) ferromagnet is determined by the spin waves which 
represent the Goldstone bosons of this nonrelativistic system. These low-energy excitations 
obey a quadratic dispersion relation. As is well-known, the spin waves of an anti ferromagnet, 
on the other hand, follow a linear dispersion law. Accordingly, the low-temperature prop- 
erties of these two systems are quite different - the difference reveals itself, e.g., in the 
low-temperature expansion for the corresponding order parameters. In this section we are 
going to address the two systems within the effective Lagrangian perspective. 

The low-temperature analysis of an 0(N) antiferromagnet, exhibiting a spontaneously 
broken internal symmetry 0(N) — > O(N-l), was the object of reference.0 Here, we only want 
to point out some basic ingredients in order to compare the system with the ferromagnet. 

In the underlying theory, the 0(N) symmetry of the Heisenberg model is explicitly broken 
by an external anisotropy field h. It is convenient to collect the (N-l) Goldstone fields U a 
in a N-dimensional vector U l = (U°, U a ) of unit length, 

U\x)U\x) = 1, (7.1) 

and to take the constant external field along the zeroth axis, h l = (h, 0, . . . , 0). The Eu- 
clidean form of the effective Lagrangian for an O(N) antiferromagnet up to order p 4 then 
reads:! 

Cff f = Ipdjrdpir - E.h'U* - e l {d IJ U i d IM U i ) 2 - e 2 (d,U l d u U 1 ) 2 

+ k 1 ^(h i U i )(d tl U%U k ) - h^(VW) 2 - h^-hW. (7.2) 

In the power counting scheme, the field U(x) counts as a quantity of order one. For the 
antiferromagnet, derivatives correspond to one power of the momentum, <9 M oc p, whereas the 
external field h counts as a quantity of order p 2 . Hence, at leading order (oc p 2 ) two coupling 
constants, T and E s , occur, at next-to-leading order (oc p A ) we have five such constants, 
ei,e 2 , h, k 2 , k 3 . 

The essential point here is the fact that a term involving a single time derivative does 
not show up in the effective Lagrangian: this topological contribution is proportional to the 
spontaneous magnetization which, in the case of an antiferromagnet, vanishes. Hence, the 
corresponding equation of motion is of second order both in space and in time, its relativistic 
structure determining the number of independent magnon states: the Fourier decomposition 
contains both positive and negative frequencies, such that a single real field suffices to 
describe one particle. Accordingly, there exist two different types of spin-wave excitations 
in an antiferromagnet - as it is the case in Lorentz-invariant theories, Goldstone fields 
and Goldstone particles are in one-to-one correspondence. These low-energy excitations 
follow a linear dispersion relation. As a consequence, in the power counting scheme for 
antiferromagnetic magnons, powers of momentum are on the same footing as powers of 
energy or temperature: k 2 oc u 2 ,T 2 . 

Let us now turn to the evaluation of the partition function of an O(N) antiferromagnet, 
N >2. The relevant graphs are depicted in Figure 4. Lorentz invariance ensures that only 
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6a 6b 6c 

FIG. 4. Feynman graphs occurring in the low-temperature expansion for the partition function 
of an O(N) antiferromagnet up to order p & . The numbers attached to the vertices refer to the 
piece of the effective Lagrangian they come from. Vertices associated with the leading term of 
£-^ff are denoted by a dot. The numbers specifying individual graphs correspond to the power of 
momentum in the derivative expansion. Note that, in the counting scheme for antiferromagnetic 
magnons, two powers of momentum correspond to two powers of energy or temperature. 

even powers of momentum occur and that loop graphs are suppressed by two powers of 
momentum. Shown are all contributions to the free energy density up to order p 6 or, 
equivalently, T 6 . Note the striking difference with respect to Figure 2, which displays the 
Feynman graphs relevant for the evaluation of the partition function of an 0(3) ferromagnet. 

The tree graphs 4b and 6c are temperature independent and hence merely renormalize 
the vacuum energy. One-loop graph 6b exclusively contains a vertex which is quadratic 
in the Goldstone boson field and thus only contributes to a renormalization of the mass 
M 2 = Yi s h/T 2 in the free Bose gas term related to graph 4a. The final expression for the 
renormalized mass takes the form0 



Mi 



Y, s h N-3 (E s hf 



JT2 



+ 



32vr 2 J^ 6 



In A + Q(h 3 ). 



h 



(7.3) 



M 



The logarithmic scale hu is determined by k± and /c 2 , i.e., by two next-to-leading order 
coupling constants (order p 4 ). The occurrence of such "chiral logarithms" is characteristic of 
the effective Lagrangian method in the Lorentz-invariant regime. To the order p 6 considered 
here, we are thus left with only one candidate for the interaction: two-loop graph 6a. 

In the limit of a zero external field, the low-temperature expansion of the pressure 
amounts to 



P 



AF 



±n 2 (N-l)T 4 + (9(T 8 



(7.4) 
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The T 4 contribution represents the free Bose gas term which originates from one-loop graph 
4a. The effective interaction among the Goldstone bosons, remarkably, only manifests itself 
through a term of order T 8 , which is beyond the scope of the present discussion. As it 
turns out, the contribution of order p 6 coming from two- loop graph 6a, is proportional to 
the external field, such that, in the limit h — *■ considered here, it vanishes. 

It is interesting to note that this two-loop graph, which involves a four magnon-vertex 
from the leading Lagrangian C^t, does not contribute to the pressure of an 0(3) ferromag- 
net, either. The reason, however, is quite different. For the antiferro magnet, the contribution 
is only absent in zero external field, whereas for the ferromagnet the contribution from graph 
8 vanishes because of space rotation symmetry or parity. 

Let us now turn to the order parameters. In the case of an 0(3) antiferromagnet this 
quantity is referred to as staggered magnetization. It is given by the logarithmic derivative of 
the partition function with respect to the external field h and, for an 0(N) antiferromagnet, 
the low-temperature expansion amounts to (for details se e§) 

K[T) I K{0) = <™^> % + ( r ) . (7.5) 

On the other hand, the low-temperature expansion for the spontaneous magnetization of 
the 0(3) ferromagnet, derived in the present paper, reads 



£(T)/£(0) = 1 - a T2 - a x T* - a 2 T* - a 3 T 4 + 0(T 



Comparing these formulae for the staggered and spontaneous magnetization, respectively, 
the following issues are worth pointing out. 

Obviously, as far as the structure of the low-temperature expansion is concerned, the 
two series are quite different. The leading terms (oc T 2 for the antiferromagnet, oc T 3 / 2 for 
the ferromagnet) are an immediate consequence of the corresponding dispersion laws: linear 
for the antiferromagnet, quadratic for the ferromagnet. They are both related to a one- loop 
graph. 

Now we have to emphasize that the evaluation of the partition function for the antiferro- 
magnet was performed in a Lorentz-invariant framework, exhibiting a relativistic dispersion 
law 

^ = v *M 2 + v 2 k\ (7.6) 

with the speed of light replaced by the spin-wave velocity v. The point is that we did 
not consider corrections proportional to k 4 ,k 6 , . . . in the dispersion relation in order to 
account for the discreteness of the lattice. Equivalently, we did not consider terms in the 
subleading pieces £ 4 ^,£^, ... of the effective Lagrangian which, although space rotation 
invariant, are no longer Lorentz invariant. They would have manifested themselves as terms 
of order T 4 , T 6 , ... in the staggered magnetization of the antiferromagnet. Like the analogous 
terms of order T 5 / 2 , T 7//2 , . . . in the spontaneous magnetization of the ferromagnet, these 
contributions merely correspond to noninteracting magnons. 

What is remarkable, however, is the fact that the two-loop contribution involving a four- 
magnon vertex from the leading piece C? e ^ of the effective Lagrangian vanishes in either 
case: in the low-temperature series for the staggered and spontaneous magnetization, the 
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spin-wave interaction does not manifest itself through this graph. Note that the reasons are 
different. For an O(N) antiferromagnet, there is actually a term proportional to T 4 in the 
low-temperature series ( |7.5| ) for the order parameter. But for the particular case N=3 we 
are considering, the term does not show up in the staggered magnetization - corrections 
associated with the magnon interaction only appear at order T 6 .EH On the other hand, the 
two- loop contribution for the ferro magnet, proportional to T 3 , vanishes due to space rota- 
tion symmetry - corrections originating from the magnon interaction manifest themselves 
through the Dyson term proportional to four powers of the temperature. 

Another striking difference concerns the fact that for the antiferromagnet, the low- 
temperature expansion is completely fixed by a single coupling constant JF up to order 
T 4 , i.e., including the first correction coming from the interaction.il This is not the case for 
the ferromagnet, where the leading interaction term, i.e., the two-loop contribution a^T 4 , 
involves coupling constants from the next-to- leading order Lagrangian £e//- We conclude 
that, for Lorentz-noninvariant systems, kinematics is much less restrictive - less information 
is available via symmetry considerations, such that more phenomenological input for the 
corresponding nonrelativistic system is required. 

VIII. CONCLUSIONS AND OUTLOOK 

The low-energy behavior of an 0(3) ferromagnet is determined by its low-energy excita- 
tions, the magnons or spin waves, which represent the Goldstone bosons of the spontaneously 
broken internal symmetry 0(3) — > 0(2). The system may be analyzed within the framework 
of nonrelativistic effective Lagrangians, which tackles the phenomenon of spontaneous sym- 
metry breaking from a unified point of view. The method exploits the symmetry properties 
of the underlying theory, i.e., the Heisenberg model in the present case, and formulates 
the dynamics in terms of Goldstone boson fields. At large wavelengths, the microscopic 
structure of condensed matter systems does not play a significant role: in the corresponding 
effective Lagrangian, the specific properties of the system only manifest themselves in the 
numerical values of a few coupling constants. 

The low-energy excitations of an 0(3) ferromagnet obey a quadratic dispersion relation. 
In the momentum power counting scheme, which provides us with a systematic expan- 
sion of physical quantities in powers of inverse wavelength, two powers of momentum thus 
correspond to only one power of energy or temperature. Remarkably, unlike in a Lorentz- 
invariant framework where loop corrections are suppressed by two powers of momentum, 
loops involving ferromagnetic magnons are suppressed by three momentum powers. 

In this nonrelativistic framework, the evaluation of the partition function for an 0(3) 
ferromagnet is presented up to order p 10 : up to order T 5 in the free energy density or, equiv- 
alently, up to order T 4 in the spontaneous magnetization. In agreement with Dyson's micro- 
scopic analysis, we find that, in the spontaneous magnetization, the spin- wave interaction 
only starts contributing at order T 4 . Moreover, within the effective Lagrangian perspective 
we readily understand the absence of a T 3 term in the spontaneous magnetization. This 
spurious term, which was the object of various microscopic investigations, vanishes because 
of space rotation symmetry or parity of the leading order effective Lagrangian. The effective 
Lagrangian method not only proves to be more efficient than the complicated microscopic 
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analysis, but also addresses the problem from a model-independent point of view based on 
the symmetries of the system. 

It is very instructive to compare the 0(3) ferromagnet and the 0(3) antiferromagnet 
within the framework of effective Lagrangians. Spin waves in antiferromagnets obey a linear 
dispersion relation, implying that powers of momentum, energy and temperature are all 
on the same footing: the power counting is identical with the standard, Lorentz-invariant 
scheme. Accordingly, the low-temperature properties of these two systems are quite different 
as illustrated by the low-temperature expansion for the pressure and for the corresponding 
order parameters, the spontaneous and staggered magnetization, respectively. 

Having established the power counting scheme for the ferromagnet, we have paved the 
way for further investigations of this nonrelativistic system within the effective Lagrangian 
perspective. In fact, the low-temperature expansion for the partition function has already 
been carried to order p 11 , where the first three-loop graphs show up. An outline of the calcu- 
lation^which goes beyond Dyson's microscopic analysis, will be presented in a forthcoming 
paper £3 
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